home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / SPPDI.z / SPPDI
Text File  |  1996-03-14  |  2KB  |  67 lines
  1.  
  2.  
  3.  
  4. SSSSPPPPPPPPDDDDIIII((((3333FFFF))))                                                            SSSSPPPPPPPPDDDDIIII((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      SPPDI   - SPPDI computes the determinant and inverse of a real symmetric
  10.      positive definite matrix using the factors computed by SPPCO or SPPFA .
  11.  
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  14.       SUBROUTINE SPPDI(AP,N,DET,JOB)
  15.  
  16. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  17.      On Entry
  18.  
  19.      AAAAPPPP REAL (N*(N+1)/2)
  20.         the output from SPPCO or SPPFA.
  21.  
  22.      NNNN INTEGER
  23.         the order of the matrix  A .
  24.  
  25.      JJJJOOOOBBBB INTEGER
  26.         = 11   both determinant and inverse.
  27.         = 01   inverse only.
  28.         = 10   determinant only.  On Return
  29.  
  30.      AAAAPPPP the upper triangular half of the inverse .
  31.         The strict lower triangle is unaltered.
  32.  
  33.      DDDDEEEETTTT REAL(2)
  34.         determinant of original matrix if requested.
  35.         Otherwise not referenced.
  36.         Determinant = DET(1) * 10.0**DET(2)
  37.         with  1.0 .LE. DET(1) .LT. 10.0
  38.         or  DET(1) .EQ. 0.0 .  Error Condition
  39.  
  40.      AAAA division by zero will occur if the input factor contains a zero on the
  41.      diagonal and the inverse is requested.  It will not occur if the
  42.      subroutines are called correctly and if SPOCO or SPOFA has set INFO .EQ.
  43.      0 .  LINPACK.  This version dated 08/14/78 .  Cleve Moler, University of
  44.      New Mexico, Argonne National Lab.  Subroutines and Functions BLAS
  45.      SAXPY,SSCAL Fortran MOD
  46.  
  47.  
  48.  
  49.  
  50.  
  51.  
  52.  
  53.  
  54.  
  55.  
  56.  
  57.  
  58.  
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.